Heritable changes in division speed accompany the diversification of single T cell fate.

Rapid clonal expansion of antigen-specific T cells is a fundamental feature of adaptive immune responses. It enables the outgrowth of an individual T cell into thousands of clonal descendants that diversify into short-lived effectors and long-lived memory cells. Clonal expansion is thought to be programmed upon priming of a single naive T cell and then executed by homogenously fast divisions of all of its descendants. However, the actual speed of cell divisions in such an emerging "T cell family" has never been measured with single-cell resolution. Here, we utilize continuous live-cell imaging in vitro to track the division speed and genealogical connections of all descendants derived from a single naive CD8+ T cell throughout up to ten divisions of activation-induced proliferation. This comprehensive mapping of T cell family trees identifies a short burst phase, in which division speed is homogenously fast and maintained independent of external cytokine availability or continued T cell receptor stimulation. Thereafter, however, division speed diversifies, and model-based computational analysis using a Bayesian inference framework for tree-structured data reveals a segregation into heritably fast- and slow-dividing branches. This diversification of division speed is preceded already during the burst phase by variable expression of the interleukin-2 receptor alpha chain. Later it is accompanied by selective expression of memory marker CD62L in slower dividing branches. Taken together, these data demonstrate that T cell clonal expansion is structured into subsequent burst and diversification phases, the latter of which coincides with specification of memory versus effector fate.

Early emergence of T central memory precursors programs clonal dominance during chronic viral infection.

Chronic cytomegalovirus (CMV) infection leads to long-term maintenance of extraordinarily large CMV-specific T cell populations. The magnitude of this so-called 'memory inflation' is thought to mainly depend on antigenic stimulation during the chronic phase of infection. However, by mapping the long-term development of CD8+ T cell families derived from single naive precursors, we find that fate decisions made during the acute phase of murine CMV infection can alter the level of memory inflation by more than 1,000-fold. Counterintuitively, a T cell family's capacity for memory inflation is not determined by its initial expansion. Instead, those rare T cell families that dominate the chronic phase of infection show an early transcriptomic signature akin to that of established T central memory cells. Accordingly, a T cell family's long-term dominance is best predicted by its early content of T central memory precursors, which later serve as a stem-cell-like source for memory inflation.

A scalable moment-closure approximation for large-scale biochemical reaction networks.

MOTIVATION: Stochastic molecular processes are a leading cause of cell-to-cell variability. Their dynamics are often described by continuous-time discrete-state Markov chains and simulated using stochastic simulation algorithms. As these stochastic simulations are computationally demanding, ordinary differential equation models for the dynamics of the statistical moments have been developed. The number of state variables of these approximating models, however, grows at least quadratically with the number of biochemical species. This limits their application to small- and medium-sized processes. RESULTS: In this article, we present a scalable moment-closure approximation (sMA) for the simulation of statistical moments of large-scale stochastic processes. The sMA exploits the structure of the biochemical reaction network to reduce the covariance matrix. We prove that sMA yields approximating models whose number of state variables depends predominantly on local properties, i.e. the average node degree of the reaction network, instead of the overall network size. The resulting complexity reduction is assessed by studying a range of medium- and large-scale biochemical reaction networks. To evaluate the approximation accuracy and the improvement in computational efficiency, we study models for JAK2/STAT5 signalling and NF κ B signalling. Our method is applicable to generic biochemical reaction networks and we provide an implementation, including an SBML interface, which renders the sMA easily accessible. AVAILABILITY AND IMPLEMENTATION: The sMA is implemented in the open-source MATLAB toolbox CERENA and is available from https://github.com/CERENADevelopers/CERENA . CONTACT: jan.hasenauer@helmholtz-muenchen.de or atefeh.kazeroonian@tum.de. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Inference for Stochastic Chemical Kinetics Using Moment Equations and System Size Expansion.

Quantitative mechanistic models are valuable tools for disentangling biochemical pathways and for achieving a comprehensive understanding of biological systems. However, to be quantitative the parameters of these models have to be estimated from experimental data. In the presence of significant stochastic fluctuations this is a challenging task as stochastic simulations are usually too time-consuming and a macroscopic description using reaction rate equations (RREs) is no longer accurate. In this manuscript, we therefore consider moment-closure approximation (MA) and the system size expansion (SSE), which approximate the statistical moments of stochastic processes and tend to be more precise than macroscopic descriptions. We introduce gradient-based parameter optimization methods and uncertainty analysis methods for MA and SSE. Efficiency and reliability of the methods are assessed using simulation examples as well as by an application to data for Epo-induced JAK/STAT signaling. The application revealed that even if merely population-average data are available, MA and SSE improve parameter identifiability in comparison to RRE. Furthermore, the simulation examples revealed that the resulting estimates are more reliable for an intermediate volume regime. In this regime the estimation error is reduced and we propose methods to determine the regime boundaries. These results illustrate that inference using MA and SSE is feasible and possesses a high sensitivity.

CERENA: ChEmical REaction Network Analyzer--A Toolbox for the Simulation and Analysis of Stochastic Chemical Kinetics.

Gene expression, signal transduction and many other cellular processes are subject to stochastic fluctuations. The analysis of these stochastic chemical kinetics is important for understanding cell-to-cell variability and its functional implications, but it is also challenging. A multitude of exact and approximate descriptions of stochastic chemical kinetics have been developed, however, tools to automatically generate the descriptions and compare their accuracy and computational efficiency are missing. In this manuscript we introduced CERENA, a toolbox for the analysis of stochastic chemical kinetics using Approximations of the Chemical Master Equation solution statistics. CERENA implements stochastic simulation algorithms and the finite state projection for microscopic descriptions of processes, the system size expansion and moment equations for meso- and macroscopic descriptions, as well as the novel conditional moment equations for a hybrid description. This unique collection of descriptions in a single toolbox facilitates the selection of appropriate modeling approaches. Unlike other software packages, the implementation of CERENA is completely general and allows, e.g., for time-dependent propensities and non-mass action kinetics. By providing SBML import, symbolic model generation and simulation using MEX-files, CERENA is user-friendly and computationally efficient. The availability of forward and adjoint sensitivity analyses allows for further studies such as parameter estimation and uncertainty analysis. The MATLAB code implementing CERENA is freely available from http://cerenadevelopers.github.io/CERENA/.